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Abstract 

In [5] we presented a method to recover the Maximum Entropy Density (MED) inferred from 
prices of call and digital options on a set of n strikes. To find the MED we need to numerically 
invert a one-dimensional function for n values and a Newton-Raphson method is suggested. In this 
note we revisit this inversion problem and show that it can be rewritten in terms of the Langevin 
function for which numerical approximations of its inverse are known. 

The approach of 5 is very similar to that of Buchen and Kelly (BK) [2] with the difference that 
BK only requires call option prices. Then, in [1] we presented another approach which uses call 
prices only and recovers the same density as BK with a few advantages, notably, numerical stability. 
[4] provides a detailed analysis of convergence and, in particular, gives various estimates of how far 
(in different senses) the iterative algorithm is from the solution. These estimates rely on a constant 
m > 0. The larger m is the better the estimates will be. A concrete value of m is suggested in [4], 
and this note provides a sharper value. 
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1 Introduction 

In [5] we presented a method to obtain the Maximum Entropy Density (MED) inferred from arbitrage 
free prices of call and digital options. More precisely, the option prices are assumed for a given maturity 
and a set of n strikes. A similar approach was firstly derived in [2] by Buchen and Kelly (BK) but they 
used call prices only. 

The inclusion of digital prices simplified the methodology by requiring the solution of n one-dimensional 
non-linear equations rather than the solution of BK's n-dimensional non-linear problem. The numer- 
ical method that we presented turned out to be more stable than BK's. Actually, Buchen and Kelly 
themselves mentioned the instability in their article. Orozco Rodriguez and Santosa [6] have also clearly 
pointed that out. 
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In general, when digital prices are prescribed, the MED obtained does not match BK unless the 
digital prices are the same as those implied by the BK density. The algorithm presented in our second 
work [3] does not assume that the digital prices are given and recovers the BK density without the 
instability issues. 

The key point to recovering the BK density is to find the right set of digital prices. This set is 
characterized by two equivalent conditions: 

• It is the set of digital prices that produces a continuous density. 

• It is the set of digital prices which implies the density with highest entropy. 

The latter condition suggests a clear way to recover the BK density by solving an optimization problem 
which turns out to be more stable than BK's original one. In addition, the n x n Hessian matrix of 
the objective function of our approach is tridiagonal in contrast to the full matrix that appears in BK's 
algorithm. 

Although our first algorithm does not generally produce the BK density, it is still needed because it 
is run at each step of our second algorithm. In particular, we still need to solve the n one-dimensional 
problems mentioned above. We turn our attention back to this problem. 

In the proof of Proposition 2.3 of [5] we have shown that, through a simple change of variables, 
solving the n equations is equivalent to inverting a single function at n different points. We noticed that 
this function has the profile of a probability distribution function (it is increasing and has limits and 

1 when the argument goes to — oo and +oo, respectively) and we wondered if an approximation of this 
function's inverse was already known. In this note we shall see that this is indeed the case. 

The function to be inverted is related to the Langevin function. This function also allows us to 
improve some results regarding the convergence analysis of our algorithm presented in [4]. 

2 Finding the MED through the Langevin Function 

Let Kq := < K\ < ... < K n < K n+ \ := oo. As in j4], for i — 0, n we define 



The function c„ has a different treatment and, as seen in [5], its derivative is easily inverted. Therefore, 
we focus on the case i < n. 

Recall that the apparent singularity at j3 = is not effective and Ci is twice continuously differentiable. 




for i < n and (3 ^ 0, 



for i < n and j3 = 




for i — n and /3 < 0, 
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For each i < n, let U := ^(K i+ i + K{) and Vi := ^(K i+ i — Ki). For i < n and /3 ^ we have 



<k(P) = In 



= ln ( e «tf«+V«) _ e/W-vO) _ lnj g 



/3J7i + In 2 



ln/3 = jSi/i +ln2 + In (sinh(^)) ~ ln /?- 



Therefore, 



= + v < Spy - £ = u * + v * ( coth ^ - w) = Ui + ViL{m > 

where L(x) := coth(a;) — a; -1 is the Langevin function. The graph of L is shown in Figure [TJ 
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Figure 1: Graph of the Langevin function. 
Recall that finding the MED, given call and digital prices, requires solving for /3j 

c'iWi) = Ki 

where Ki G (Ki, JQ+i) is given in terms of digital and call prices only. Hence, we need to solve 

Ui + V t L(pV t ) = K 

which gives 

_ 1_ ( Ki - Uj 

Pl ~Vi L \ Vi 

A few approximation formulas for L _1 have been proposed in the literature. Firstly, through the 
Taylor expansion around x = 0, 

r . 9 3 297 5 1539 7 

L 1 (x) = 3a; + -ar H x b H x 1 + ... 

V ; 5 175 875 
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Notice that the domain of L _1 is (— 1, 1) and L (x) — > ±00 when x — > ±1. Therefore a growing number 
of terms of the Taylor expansion is needed when |x| — > 1. To overcome this issue, from the Taylor 
expansion above one can derive the following Pade approximant [5] 

_ 36„,2 

L-\x)=x 6 -^t- 2 +0{x% 

1 - m x 



However, the poles of this approximation are located at x = ±y35/33 and not at ±1 as those of L 1 . 
Finally, rounding 36/35 « 1 and 33/35 « 1 yields the rounded Pade approximant [3] 

L- 1 (x)=x-' 



1 -x 2 

Another approximation formula with a relative error smaller than 6.4 • 10 -4 [I] is given by 

_ J 1.31446 -tan(l. 58986a;) + 0.91209a;, if \x\ < 0.84136, 

^ ~~ I l/(sign(a;) - a;), if 0.84136 < |x| < I. 



3 Convergence Analysis 

We recall some notation from [3] before stating stronger results for the convergence analysis presented 
in Section 6.2 of that article. 

For a given maturity, we fix a set of undiscounted arbitrage- free call prices Ci, C n at prescribed 
strikes K\, K n and also the forward asset price Co (which is seen as the price of a call option with 
strike Kq :— 0). These prices impose a number of constraints on the values that undiscounted arbitrage- 
free digital prices D\, D n on the same strikes can assume. Let ft be the set of D := (Di, D n ) that 
verify these constraints. 

For each D G 57 the results of [5 show that there exists a density which maximises entropy 
subject to implying the same prices of calls Co, C\, C n and digitals D\, D n . Let H{D) be the 
entropy of such a density. By varying D in order to maximise H, we recover the BK density and its 
corresponding vector of digital prices D. 

Proposition 6.1 of [3] introduces a constant m > that drives the subsequent convergence analysis. 
More precisely, it is shown that for all D e Q we have 

(i) . H{b)-H{b)<^\\H'{D)f; 

(ii) . ||D-D||<£||Jf'(£)||;and 
(hi). \\ 9b - 9£) \\ L1 <-L=\\H'(b)\\. 

From the relations above, we see that the bigger m is, the better our knowledge of the convergence 
becomes. In other words, the number of steps required to get an approximated solution within a given 
tolerance decreases. Proposition 6.1 also suggests taking m = 4 sin 2 (7r/ (2n+2)), but this value approaches 
as the number of strikes goes to infinity. In the following, we shall see that we can improve this result 
and take m = 4sin 2 (7r/(2n + 2)) + 1/2. 
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Looking at the proof of Proposition 6.1, we see that m is given by the sum of m\ := 4sin 2 (7r/(2n + 2)) 
and 

= 1 . f (gi - Xp) 2 . f (gj - i^) 2 QKj+i - if t ) 2 ] (K n - K n f 

1712 ■ 2 mm \ p oC [,'(A>) ' i=W-i\ PicJ'dSj) ' ftcJ'GSi) /' p„<(/3 n ) 

where pi = Di — ZV|_i and X; = c4(/3i) for all i = 0, ...,n. We shall see that 7712 > 1/2 and for this, it 
suffices to show that 



fa) 

and 



> 1 Vi = 1, n, 



(iT z+1 -c^)) 2 

S^j - 1 Vj = °»-' B - L 

The arguments are similar and we shall show only the first relation above. 
Firstly, we notice that pt < 1; therefore, it is enough to show that 

(4(ft) Vi=1 „ 

dim - '-' 

The case i = n is straightforward and is left to the reader. In the sequel we assume i < n. 
Using (JTJ and Ui - Ki = Vi yields 

(cJCSj) - i^) 2 = (C/ 4 + ViL(piVi) - X,) 2 - (1 + ^GW) 2 

and 
Hence 

( C ;(ft)-^) 2 = (l + Ljm)) 2 

To finish the proof we shall show that 

{1 + L(x)) 2 >L'(x) VxER. 
For x = the result follows from L(0) = and 1/(0) = 1/3. For a; 7^ 0, the result is equivalent to 

(l + coth a; — a; 1 ) 2 > 1 — coth 2 x + x~ 2 coth 2 x + 2 coth a; — 2x~ 1 — 2a?~ 1 coth a; > — coth 2 : 

<==> coth 2 x + coth a; — x^ 1 — x^ 1 coth a; > 
(coth a; - aT 1 ) (coth a; + 1) > 
L{x) (cothx + 1) > 0. 

The last inequality follows from the fact that L(x) and coth a; + 1 both have the same sign as x. 
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